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At the precision reached in current lattice QCD calculations, electromagnetic effects are becoming 
numerically relevant. Here, electromagnetic effects are included by superimposing U(l) degrees 
of freedom on Nf = 2 + 1 QCD configurations from the Budapest-Marseille-Wuppertal Collabo- 
ration. We present preliminary results for the electromagnetic corrections to light pseudoscalars 
f — mesons masses and discuss some of the associated systematic errors. 

(N 

o 

(N 



X 



The XXIX International Symposium on Lattice Field Theory - Lattice 2011 
July 10-16, 2011 

Squaw Valley, Lake Tahoe, California 



* Speaker. 

+ Unite Mixte de Recherche (UMR 7332) du CNRS, de l'Universite d' Aix-Marseille et de 1'Universite Sud Toulon 
Var. Unite affiliee a la FRUMAM, Federation de Recherche 2291 du CNRS. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http : / /pos . sissa . it / 



Systematic errors in partially-quenched QCD plus QED lattice simulations 



A. Portelli 



1. Motivation 



Isospin is a near symmetry of the hadron spectrum because it is only broken by small effects: 

(i) the mass difference m u — m^ 

(ii) the difference in the charge of the u and the d quark 
which are summarized in the following table 1 : 





u 


d 


Mass (MeV) [1,2] 
Charge 


2.15(03) stat .(10) sys . 

h 


4.79(07) stat .(12) sy , 



These effects are expected to be at the percent level. The size of mass breaking is the mass differ- 
ence m u — m c \ relatively to a typical QCD scale Aqcd and the order of electromagnetic breaking 



is the fine structure constant at zero momentum a 



An 



137- 



However, their importance is not 



commensurate to their size. For instance, they are responsible for the stability of matter through 
the proton-neutron mass difference. Moreover, lattice QCD calculations have recently approached 
percent or even sub-percent precision [4, 5], so the inclusion of these effects becomes relevant. 
Another interesting isospin breaking quantity is the absolute correction to Dashen 's theorem 



where 



AaD 



a em m| 



AemM|-A E mM^ 



(M| 4 



-M; 



pO)m„=m d 



(1.1) 



(1.2) 



is the electromagnetic squared mass splitting of the isospin multiplet P. One can also consider the 
dimensionless relative correction to Dashen 's theorem : 



A E mM2 



(1.3) 



R. Dashen has shown in [6] than AaD = in the SU(3) chiral limit and that the leading corrections 
are 0(am s , a 2 ). The quantity AaD is interesting because it is very sensitive to the up and down 
quark masses. Moreover, a precise estimation of these quantities has not been made until now 
(cf. Table 1). 



2. Simulation setup 

For this work, we used a subset of our 2010 SU(3) gauge configurations [2]. These configura- 
tions were generated using Nf = 2 + 1 QCD simulations with the tree-level Symanzik gauge action, 
tree level 0(a) -improved Wilson fermions and two steps of HEX smearing. We have already used 
this dataset to compute light quark masses [1, 2] and the kaon bag parameter [18]. 

'We have chosen to quote the quark masses from [1,2]. Please see [3] for a complete list of lattice determinations 
of these masses. 
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A A D (MeV z ) 


A R D 






1230 


0.80 


Donoghue et al. [7] (1993) 


00 
O 

"o 


1300(400) 


1.0(3) 


Bynens [8] (1993) 


<u 


360 


0.26 


Baur and Urech [9] (1996) 


s 

o 
ss 


1060(320) 


0.9(4) 


Bynens and Prades [10] (1997) 


<D 


1080 


0.68 


Gao et al. [11] (1997) 




1070 


0.74 


Bynens and Damelsson [12] (2007) 




526 


0.39 


Duncan et al. [13] (1996) 




340(92) 


0.30(8) 


Blum et al. [14] (2007) 




1250(550) 


N/A 


Basak et al. [15] (2008) 




830(180) 


0.60(14) 


Portelli et al. [16] (2010) 




707(75) 


0.63(6) 


Blum et al. [17] (2010) 



Table 1: Results for the violations to Dashen's theorem in phenomenelogy and lattice computations. Order 
is chronological. Bold numbers are the results given by the authors, the others are deduced from information 
given in the corresponding paper. 

To include QED, we generate real electromagnetic fields using a non-compact formulation 
(cf. [16] for more details). Then we phase SU(3) strong links by the U(l) links exp(;'2A |U ), where 
Q is a chosen electric charge. The resulting U(3) links are used inside the Dirac- Wilson operator 
to compute quark propagators. This method allows to use previously generated SU(3) gauge con- 
figurations to obtain results, but the electromagnetic vacuum polarization is not taken into account. 
For the preliminary study presented here, the masses of the up and down valence quarks are tuned 
individually such that the bare masses are equal. This occurs when these masses become equal to 
the light sea quark mass [16]. The n° squared mass is obtained by averaging squared ground state 
energies obtained with uu and dd connected pseudoscalar correlators. This is correct up to NLO 
isospin breaking corrections. 

In this setup, one is confronted with two new types of systematics effects : the quenching of 
the electromagnetic field and electromagnetic finite volume effects. 

3. QED quenching errors 

Using photon fields without vaccuum polarization leads to partial quenching effects: valence 
quarks have electric charges but sea quarks do not. These effects can be evaluated in partially 
quenched chiral perturbation theory with QED (PQ^PT+QED). Next to leading order (NLO) 
SU(3) formulas for self-energies and decay constants can be found in [12, 17]. The SU(3) NLO 
sea contribution to a pseudoscalar meson squared mass is given by : 

-wYiXijt^ (3-D 

3 s=4 

where i and j are the valence indices of the quarks composing the meson, indices, s, between 4 
and 6 are sea quark indices, Xkl = Bo{mk + mi), e ~ 0.302822 is the positron electric charge, q\ are 
quark charges in units of e and Y\, C and Fq are low energy constants (LECs). 



e l C 
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For the results presented here, the sea and valence masses are very nearly equal and we can 
assume that the absence of sea charges is the only partial quenching effect. In that case, one can 
check easily that (3.1) is independant of the scale pL. The LECs C and Fq are essentially known [12], 
but Y\, which is a sea electromagnetic contribution, is unknown. RBC, who studied other partially- 
quenched QED LECs [14, 17], considers that it would be unnatural that Y\ > 10~ 2 . Thus, in the 
following, we will use Fi = 1(T 2 . With (3.1), one obtains the following partial quenching error 
estimates: 

• negligible (less than 0.1%) for M K + , M K o , A E m^ and A EM M 2 

• ~0.4%forM 3r o 

• ~ 5% for A EM M/s: and A EM M^ 

4. Electromagnetic finite volume effects 

It is already known [2, p. 18] that the Budapest-Marseille-Wuppertal gauge ensembles have 
spatial volumes that allow to neglect QCD finite volume effects. However, as electromagnetism is 
a long range interaction, one might expect large QED finite volume effects on a periodic lattice. 
Predictions for these effects have been made in PQ^PT+QED [19, 17]. These formulas are far 
from simple and complicate chiral fits. Additionally, their predictivity seems limited [17, p. 16]. 

Here we choose a more straightforward approach. Dimensional analysis suggests that the 
leading finite volume correction to the splitting of a squared hadron mass, AM 2 , takes the form : 



where L is the spatial extent of the lattice and A is an unknown dimensionless constant. We found, 
as presented in the next section, that finite volume corrections of type (4.1) fit lattice data well but 
are not compatible with SU(3) PQ^PT+QED (cf. Figure 1). 

5. Preliminary results 

For this preliminary analysis, we used the /3 = 3.31 subset of Budapest-Marseille-Wuppertal 
2-HEX gauge configurations [2]. The main features of this dataset are: 

• one lattice spacing a ~ 0. 12 fm 

• 12 pion masses from 135 MeV to 422 MeV 

• 4 spatial volumes from (2 fm) 3 to (5.8 fm) , with M n oL > 4 

To interpolate a quantity to the physical masses, we use a Taylor expansion in M 2 and 
M\ = ^(M^ + +M^ — M 2 + ). We choose M 2 because it has negligible finite volume corrections 
compared to those in M 2 + . The lattice spacing is computed in physical units using the Q. baryon 
mass. The infinite volume extrapolation is performed by including corrections of the kind (4. 1) into 
the fit. Fits are earned out using fully correlated % 2 minimization and a bootstrap error analysis. 



AM 2 H-AM 2 (L)=e : 



2 A 
L 2 



(4.1) 
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Finite volume effects on AM% 
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Figure 1: Comparison of our results with different models of finite volume corrections for the kaon squared 
mass splitting, as a function of lattice spatial extent. The red band represents the fit of our lattice data to 
the model (4.1), with M K o fixed to its physical value. Only statistical errors are shown. The green curve 
represents the prediction of SU(3) PQ^PT+QED, which depends only on the LO, LECs Fo and C, given 
in [12]. There is a clear disagreement between the two. However, the jj contribution from the full SU(3) 
PQjPT+QED expression (blue curve) roughly agrees with our data. 



Using this methodology, we obtain the following preliminary results (an example of such a fit is 
shown in Figure 2) : 

A EM M 2 K = 2179(34) stat .(100) qu .(?) sys . MeV 2 (5.1) 
A EM M 2 = 1283(30) stat .(0) qu .(?) sys .MeV 2 (5.2) 
A A D = 896(37) stat . (100) qu . (?) sy , MeV 2 (5.3) 
A R D = 0.70(4) sta ,(8) qu .(?) sys . (5.4) 

where (?) sys . stands for the systematic errors which we have not yet estimated, such as those asso- 
ciated with taking the continuum limit, etc. 



6. Conclusion 

Using the methodology that we presented last year [16], we extended our dataset down to the 
physical value of the light quark mass. We also studied two important systematics effects: QED 
quenching and finite volume effects. 

Our results are promising and close to phenomenological estimates [3]. Moreover, concerning 
the finite volume effects, the results are compatible with a simple jj model and it seems that SU(3) 
PQ^PT+QED fails to describe them. 
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In the short term, we will continue our analysis on several lattice spacings with non-degenerate 
up and down quark masses. This will provide the last ingredients to define precisely the physical 
point, which require a continuum limit and physical isospin breaking. 
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Figure 2: Pion squared mass splitting vs. neutral pion squared mass. Red points represent raw lattice data, 
dark red points are the same data with finite volume correction of type (4.1) and the dashed blue line is the 
result of the physical point fit. 
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